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ABSTRACT 

We compare the non-linear matter power spectrum in real space calculated analytically from 3rd- 
order perturbation theory with iV-body simulations at 1 < z < 6. We find that the perturbation theory 
prediction agrees with the simulations to better than 1% accuracy in the weakly non- linear regime 
where the dimensionless power spectrum, A^(fc) = fc^P(fc)/27r^, which approximately gives variance 
of matter density field at a given fc, is less than 0.4. While the baryonic acoustic oscillation features 
are preserved in the weakly non-linear regime at z > 1, the shape of oscillations is distorted from the 
linear theory prediction. Nevertheless, our results suggest that one can correct the distortion caused 
by non-linearity almost exactly. We also find that perturbation theory, which does not contain any 
free parameters, provides a significantly better fit to the simulations than the conventional approaches 
based on empirical fitting functions to simulations. The future work would include perturbation theory 
calculations of non-linearity in redshift space distortion and halo biasing in the weakly non-linear 
regime 



Subject headings: cosmology : theory — large-scale structure of universe 



1. INTRODUCTION 

Cosmological linear perturbation theory has been re- 
markably successful in explaining the precision mea- 
surements of temperature and polarization anisotropies 
of the cosmic microwave background (CMB), most no- 
tably from the Wilkinson Micr owave Anisotropy Probe 
(WMAP) IjBennett et al. 2003|) . The CMB data, com- 
bined with linear theory, have enabled us to determine 
many of t he cosmological para meters to better than 10% 
accuracy (|SDergel et al. 2006|) . As the CMB data im- 
prove, however, it has become increasingly clear that one 
has to combine the CMB data with the other probes to 
break degeneracies between the parameters that cannot 
be constrained very well by the CMB data alone. For ex- 
ample, the CMB alone cannot break degeneracy between 
the equat ion of state of dark energy, w, and matter den- 
sity, n^ (jSpergelet al. 2006J) . 

The large-scale structure (LSS) of the universe has 
also been known as an excellent probe of cosmologi- 
cal fluctuations as well as cosmological parameters, as 
proven successfully by the Two-decree Fiel d Galaxy 
Redshift Survey {2dFGRS) llCole et al. 20051) and the 
Sloan Disital Skv Survey (SPSS) llTegmaTk et ah 2001 
iSeliaket al. 2005J) . A joint analysis of the future CMB 
and LSS data is extremely powerful in constraining most 
of the cosmological parameters to better than a few per- 
cent accuracy (e.s.. iTakada et al. 200(Tt) . In particular, 
the LSS data would allow us to constrain "additional" 
parameters such as the mass of neutrinos and the shape 
of the primordial power spectrum, which would remain 
relatively poorly constrained by the CMB data alone. 

The success of this approach depends on our ability to 
predict the power spectrum of CMB and LSS from the- 
ory. Linear theory provides adequate precision for CMB, 
as the amplitude of CMB anisotropy is only 10~^; how- 
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ever, theory of LSS has not reached to the point where 
one can use LSS for precision cosmology at the level sim- 
ilar to CMB. There is a larger degree of non-linearity 
in LSS. In order for the LSS data to be as powerful as 
the CMB data, it is crucial that we can predict the LSS 
power spectrum to 1% accuracy, which is nearly one or- 
der of magnitude better than the current precision. 

In principle, theory of LSS may be developed using A^- 
body simulations. This approach has been widely used 
in the literature. One method builds on the so-called 
HKLM formalism l)Hamilton et al. 19911) , which interpo- 
lates between the linear regime on large scales and the 
stable clustering regime on small scales using a fitting 
function to A^-body simulatio ns. The HKLM me thod 
was further elaborated by IjPeacock fc Doddsl996^1 . 
The other method builds on the so-called halo model 
IjScherrer fc Ber tschingcr 1991*1. which w as furth er elab- 
orated by e.g., (Seliak 2000; Smith ct al . 20031) . Both 
approaches are based on empirical methods, fitting to 
A^-body simulations mainly at z ~ 0. While these predic- 
tions may be good to within 10%, one should not expect 
1% accuracy from these. Also, these methods, in their 
current form, do not allow for non-linearity in redshift 
space distortion in the weakly non-linear regime, which 
limits their practical use for the actual data analysis. 

We use an alternative approach based on cosmological 
perturbation theory (PT). One can calculate the next-to- 
leading order correc tion to the linear power spectrum by 
using 3rd-order PT IjVishniac 1983": Suto & Sasaki 1991; 
Makino et al. 19921: Jain fc Bertschingcr 1994 : 



Scoccimarro fc Frieman 1996|) . The advantage of 

PT is that it provides an exact solution for the non- 
linear matter power spectrum as long as one applies it 
to the region in k space where perturbative expansion is 
valid. (We shall call this region in k space the "weakly 
non- linear regime".) One still needs to use simulations 
to find the maximum k below which perturbation 
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expansion is valid, which is one of the goals of this 
paper. 

Cosmological PT, including non-linear corrections to 
the power spectrum, wa s actively investigated in 1990's 
IjBcrnardc au et al. 2002. for a review). In particular, a 
lot of efforts have been devoted into understanding the 
non-linear power spectrum at z ~ 0. It was shown that 
perturbation approach would not provide accurate de- 
scriptions of the power spectrum at z ~ due to too 
strong non-linearity. Our results are consistent with the 
previous work; however, we focus on the power spectrum 
at z > 1, where non-linearity is still modest and thus PT 
should perform better. 

Our work is motivated by recent proposals of high- 
z galaxy survey proje cts such as the Cosmic Inflation 
Probe (CIP) ( Melnick eTal. 2005fl. Hobby-Ebery Dark 
Energy Experiment (HETDEX) l|Hill et al. 2004|) . and 
Wide-Geld Fi ber-fed Multi Object Sp ectrograph {WF- 
MOS) survey IjGlazebrook et al. 2005J) . to mention a few. 
The goal of these missions is to measure the power spec- 
trum of high-z galaxies to a few percent accuracy. These 
missions should be able to measure the baryonic fea- 
tures in the power spectrum accurately. On the other 
hand, it has been pointed out that non-linearity would 
distort the baryonic features in a complex way so that 
it might be challenging to extract the underlying bary- 
onic features from th e observed galaxy power spectrurn 
lIMeiksin et al. 19991: Is^ringel et al. 2005t IWhite 200^ 
iSeo fc Eisenstein 200,'t() . We show that, as far as non- 
linearity in the matter power spectrum in real space is 
concerned, we can correct it almost exactly. 

This paper is organized as follows. We briefly review 
the 3rd-order PT in § 2, and describe our A^-body nu- 
merical simulations in § 3. We compare the analytical 
predictions with simulations in § 4. We pay a particular 
attention to non-linearity in the baryonic acoustic oscil- 
lations. We give discussion and conclusions in § 5. We 
test convergence of our results in Appendix ^ 

2. NON-LINEAR MATTER POWER SPECTRUM: 
3RD-ORDER PERTURBATION THEORY 

We review 3rd-order PT calculations of the 
next-to-leading order correction to the mat- 
ter power spectrum, following the pioneer- 
ing wo rk in the liter atu re d Vishniac 1983^ 
Frv 19841: IGoroff et al.T98^ ISuto fc Sasaki 1991; 
Makino et al. 1992 | iJain fc Bertschinger 19941 

Scoccimarro fc Frieman 199(tD . As the power spec- 
trum, P{k,T), is a quadratic quantity of the density 
field in Fourier space, ^k ('''), 

(4(r)5l(r)) = {27TfP{k,T)Sn{k--k'), (1) 

the 3rd-order expansion in the density field is neces- 
sary for obtaining the next-to-leading order correction 
to P(fc,r). We often use the "dimensionless power spec- 
trum" , A^ (fc, r), which represents the contribution to the 
variance of density field per In fc. 



(<52(x,r)) = 



dk 



A2(fc,r), 



(2) 



where A^{k,T) = k^P{k,T)/{2n'^). 

We treat dark matter and baryons as pressureless dust 
particles, as we are interested in the scales much larger 



than the Jeans length. We also assume that peculiar ve- 
locity is much smaller than the speed of light, which is 
always an excellent approximation, and that the fluctu- 
ations we are interested in are deep inside the horizon; 
thus, we treat the system as Newtonian. The basic equa- 
tions to solve are given by 



6 + V-[{l + d)v] =0, 



V + (v • V)v = --V - V(/), 
a 



(3) 



(4) 



V^(j> = AnGa^pS, (5) 

where the dots denote d/dr (t is the conformal time) , V 
denotes d/dx (x is the comoving coordinate), v = dx/dr 
is the peculiar velocity field, and cj) is the peculiar grav- 
itational potential field from density fluctuations. We 
assume that v is curl-free, which motivates our using 
9 = V • V, the velocity divergence field. Using equa- 
tion |SJl and the Friedmann equation, we write the conti- 
nuity equation [Eq. (PJ] and the Euler equation [Eq. Q] 
in Fourier space as 



SkiT)+0k{T) 



d^ki I 3 



d-^fc2^Z5(kl +k2 -k)^i-^4,(T)^ki(r), (6) 



a 2a^ 



d^ki 

'{2^ 



d^k25D{^l + k2 - k)^^^J^-J^4i(T)0k.(T) 



ZiK-i rvn 



(7) 



respectively. 

To proceed further, we assume that the universe is mat- 
ter dominated, rim(7") = 1 and ai^r) oc r^. Of course, 
this assumption cannot be fully justified, as dark en- 
ergy dominates the universe at low z. Nevertheless, it 
has been shown that the next-to-leading order correc- 
tion to P{k) is extremely insensitive to the underlying 
cosmology, if one uses th e correct growth factor for S-^It) 
llBernarde au et al. 2002(1 . Moreover, as we are primar- 
ily interested in z > 1, where the universe is still matter 
dominated, accuracy of our approximation is even better. 
(We quantify the error due to this approximation below.) 
To solve these coupled equations, we shall expand (5k (t) 
and 9i^(t) perturbatively using the n-th power of linear 
solution, (5i(k), as a basis: 

x/j \ S^ n/ \ f d,^1i d^qn-i 



/n 
d^qnSoC^qt-'k) 



x-F'«(qi,q2,--- ,q«)'5i(qi)---'5i(q„), (8) 

Zn \ S~^ -I \ n~\i \ f d qi d g„_i 

n— 1 ^ ^ ^ ^ 



/n 
d^qn6DC^(\,-\<.) 



xG„(qi, q2, • • ■ , q„)'5i(qi) ■ ■ • (5i(q„ 



(9) 
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Here, the functions F and G as well as their r ecursion 
relations are given in l).Tain fc Bertschinger 1993) . As the 
linear density field, (5i, is a Gaussian random field, the 
ensemble average of odd powers of i5i vanishes. There- 
fore, the next-to-leading order correction to P{k) is 

P{k,T) = a^{T)Pii{k) + a\T)[2Pi:i{k) + P22{k)l (10) 
where 



P22{k) = 2 



F^'\^iM 



_ 5 



-Pii(g)Fii(|k-q|) Fi^)(q,k-q 



(11) 



2 (ki • k 
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ki ka 






kl 



fcj 



(12) 



2^i3(fc) = Sr^iiW 



252 



158- 



dq 

'9 



^11(9) 



g4 
19 42— 



+ T^^(9'-fc')^(2fc' + 7(z2)ln 



\k^q\ 



(13) 



p(s) 



While F2 (ki,k2) should be modified for different cos- 
mological models, the difference vanishes when ki || k2. 
The biggest correction comes from the configurations 
with kl ± k2, for which [i^2^"^(ACDM)/i^2^"^(EdS)]2 ~ 



1.006 and < 1.001 at z 



and z 



2 
> 



1, respec- 



tively. Here, i^2 (EdS) is given by equation (|12|l . while 

F^ ^(AC DM) contains corrections due to Vim # 1 and 
r^A ^ l|Matsubara 19951 IScoccimarro et al. 1998J) . and 
we used Vt^ = 0.27 and I^a = 0.73 at present. The in- 
formation about different background cosmology is thus 
almost entirely encoded in the linear growth factor. We 
extend the results obtained above to arbitrary cosmo- 
logical models by simply replacing a{T) in equation H10|l 
with an appropriate linear growth factor, D(z), 

P{k,z)^D^z)Pi,(k)+D\z)[2P,3{k)+P22{k)]. (14) 

We shall use equation (fTT)) - ((n|l to compute P{k,z). 

3. Af-BODY SIMULATIONS AND ANALYSIS METHOD 



We use the TVD llRvu et al. 1993|) code to simu- 
late the evolution of (^(x, r). The TVD code uses 
the Particle-Mesh scheme for gravity, and the Total- 
Variation-Diminishing (TVD) scheme for hydrodynam- 
ics, although we do not use hydrodynamics in our cal- 
culations. To increase the dynamic range of the de- 
rived power spectrum and check for convergence of the 
results, we use four box sizes, Lbox = 512, 256, 128, 
and 64 h~^ Mpc, with the same number of particles, 
N = 256^. (We use 512^ meshes for doing FFT.) We 
use the following cosmological parameters: ftm = 0-27, 
r^b = 0.043, flA = 0.73, h = 0.7, ag = 0.8, and n^ = 1. 
We output the simulation data at 2; = 6, 5, 4, 3, 2 and 1 
for 512, 256 and 128 h"^ Mpc, while only at z = 6, 5, 4 
and 3 for 64 h'^ Mpc. 

We suppress sampling variance of the estimated P{k, z) 
by averaging P{k, z) from 60, 60, 20, and 15 indepen- 
dent realizations of 512, 256, 128, and 64 h^^ Mpc sim- 
ulations, respectively. We calculate the density field on 




0.10 1.00 

Wovenumber [h Mpc"^] 

Fig. 1. — Power spectrum at 2 = 1, 2, 3, 4, 5 and 6 (from top 
to bottom), derived from A^-body simulations (dashed lines), per- 
turbation theory (solid lines), and linear theory (dot-dashed lines). 
We plot the simulation data from 512, 256, 128, and 64 h''^ Mpc 
simulations at fc < 0.24 h Mpc'^, 0.24 < k < 0.5 h Mpc"!, 
0.5 < A: < 1.4 h Mpc~l, and 1.4 < A: < 5 h Mpc"^, respectively. 
Note that we did not run 64 h~^ Mpc simulations at 2 = 1 or 2. 



512^ mesh points from the particle distribution by the 
Cloud-In-Cell (CIC) mass distribution scheme. We then 
Fourier transform the density field and average |^k(''')P 
within k — Ak/2 < |k| < /c -I- Afc/2 over the angle to 
estimate P{k,z). Here, Afc — 27r/Lbox. Finally, we cor- 
rect the estimated P{k) for loss of power due to the CIC 
pixelization effect using the window function calculated 
from 100 realizations of random p article distribution s. 

We use the COSMICS package IjBertschineer 1995J) to 
calculate the linear transfer function (with linger) and 
generate the input linear matter power spectrum and 
initial conditions (with grafic). We have increased the 
number of sampling points for the transfer function in 
k space from the default value of COSMICS, as the de- 
fault sampling rate is too low to sample the baryonic 
acoustic oscillations accurately. (The default rate re- 
sulted in an artificial numerical smoothing of the oscil- 
lations.) We locate initial particles on the regular grid 
(i.e., we do not randomize the initial particle distribu- 
tion), and give each particle the initial velocity field us- 
ing the Zel'dovich approximation. This procedure sup- 
presses shot noise in the derived power spectrum, which 
arises from randomness of particle distribution. We have 
checked this by comparing P(fc, z) from the initial con- 
dition to the input linear spectrum. However, some shot 
noise would arise as density fluctuations grow over time. 
While it is difficult to calculate the magnitude of shot 
noise from the structure formation, we estimate it by 
comparing P{k, z) from large-box simulations with that 
from small-box simulations. We do not find any evi- 
dence for shot noise at z > 1; thus, we do not sub- 
tract shot noise from the estimated P{k,z). To be 
conservative, we use 512, 256, 128, and 64 ft,^^ Mpc 
simulations to obtain P(fc, z) at fc < 0.24 h Mpc^ , 
0.24 < fc < 0.5 /i Mpc^\ 0.5 < yfc < 1.4 h Mpc~\ 
and 1.4 < k < 5 h Mpc^^, respectively, to avoid the 
residual CIC pixelization effect and potential contami- 
nations from unaccounted shot noise terms as well as ar- 
tificial "transients" from initial conditions generated by 



JEONG & KOMATSU 













. z = 5 




. 


1.0 


: Sr. 

" : Lir 


ith et ol. (2003) 

ear Spectrum^ ^^ 


::SS:S^ 


0.5 




^--"-^^'^'''^'^^^^ 


.: 


nn 


.— -— -'"'^^ 








0.0 0.2 0.4 0.6 0.8 1.0 1.2 1 .' 
Wovenumber [h Mpc"'] 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.' 
Wovenumber [h Mpc"'] 





Wovenumber [h Mpc '] 






Wovenumber [h Mpc '] 



0-0 0.2 0.4 0.6 0.8 1.0 1.2 1 .' 
Wovenumber [h Upc"'] 



Fig. 2. — (Top) Dimensionlcss power spectrum, A^(fc). The 
solid and dashed lines show perturbation theory calculations and 
A'^-body simulations, respecti vely. The dotted l ines show the pre- 
dictions from halo approach ISmith et al. 20 03). The dot-dashed 
lines show the linear power spectrum. (Bottom) Residuals. The 
errorbars show the TV-body data divided by the perturbation the- 
ory predictions minus one, while the solid c urves show the halo 
model calculations given in ISmith et al. 200a) divided by the per- 
turbation theory predictions minus one. The perturbation theory 
predictions agree with simulations to better than 1% accuracy for 
A2{fc) <0.4. 



the Zel'dovich approximation IjCrocce et al. 200a) . The 
initial redshifts are ^initial = 27, 34, 42, and 50 for 512, 
256, 128, and 64 h~^ Mpc simulations, respectively. In 
App end ix \K\ we show more on the convergence test (see 

Fig. mi). 

4. RESULTS 

Figure n compares P{k,z) at 2; = 1, 2, 3, 4, 5 and 6 
(from top to bottom) from simulations (dashed Hnes), FT 
(solid lines), and linear theory (dot-dashed lines). The 
FT predictions agree with simulations so well that it is 
actually difficult to see the difference between FT and 
simulations in Figure^ The simulations are significantly 
above the linear theory predictions at high k. 

To facilitate the comparison better, we show A^(fc, z) 
[Eq. ©l in Figure |21 We find that the FT predictions 
(thin solid lines) agree with simulations (thick solid lines) 
to better than 1% accuracy for A^(fc,z) < 0.4. On the 
other hand, the latest predictions from halo approach 




Fig. 3. — Non-linearity in baryonic acoustic oscillations. All 
of the power spectra have been divided by a smooth power 
spectrum without baryonic oscillations from equation (29) of 
JEisenstcin & Hu 1998). The errorbars show Af-body simulations, 
while the solid lines show perturbation theory calculations. The 
dot-dashed lines show the linear theory predictions. Perturbation 
theory describes non-linear distortion on baryonic oscillations very 
accurately at 2 > 1. Note that different redshift bins are not inde- 
pendent, as they have grown from the same initial conditions. The 
Af-body data at fc < 0.24 and k > 0.24 h Mpc~^ are from 512 and 
256 h~^ Mpc box simulations, respectively. 



IjSmith et a l. 2003) (dotted lines) perform significantly 
worse then FT. This result suggests that one must use FT 
to model non-linearity in the weakly non-linear regime. 
The baryonic features in the matter power spectrum 
provide a powerful tool to constrain the equation of state 
of dark energy. This method uses the fact that the CMB 
angular power spectrum sets the physical acoustic scale, 
and thus the features in the matter power spectrum seen 
on the sky and in redshift space may be used as the stan- 
dard ruler, giving us the angular diameter distance out to 
the gala xy distribution at a given survey redshift as well 
as H(z) llMa,tsuba,ra fc Sza.la'v 2003HHu fc Ha.ima.n 200,'^ 
iSeo fc Eisenstein 20031: iBlake fc Glazebrook 2003(1 . In 
order for this method to be viable, however, it is cru- 
cial to understand distortion on the baryonic acoustic 
oscillations caused by non-linearity. This has been inves- 
tigated so far mostly using direct numeric al simulations 
JMeiksin et al. Iflflflt ISnnngel et ^OOO^ IWhite 200,^ 



I 



ISeo fc Eisenstein 20 051. fMeiksi n et al. 19991) also com- 
pared the FT prediction with their A'^-body simulations 
at z = 0, finding that FT was a poor fit. This is be- 
cause non-linearity at z = is too strong to model by 
FT. Figure O shows that FT provides an accurate an- 
alytical account of non- linear distortion at z > 1: even 
at z = 1, the third peak at fc ~ 0.18 h Mpc""'^ is mod- 
eled at a few percent level. At z > 2, all the oscilla- 
tory features are modeled to better than 1% accuracy. 
A slight deficit in power from A^-body simulations at 
fc '^ 0.2 h Mpc~^ relative to the perturbation theory 
predictions at z = 2 may be due to artificial transient 
modes from the Zel'dovich approximation used to gen- 
erate initial conditions. One may eliminate such an ef- 
fect by either using a smaller box-size or a better initial 
condition from the second-order Lagrangian perturba- 
tion theory (Croccc et aL_2006). As the power spectrum 
at fc > 0.24 h Mpc~^ from 256 h~^ Mpc simulations at 
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Fig. 4. — Non-linearity and the amplitude of matter fluctuations, 
ag. In each panel the lines show the linear spectrum and non-linear 
spectrum with as = 0.7, 0.8, 0.9 and 1.0 from bottom to top. 

z = 2 agrees with the perturbation theory predictions 
very well, we conclude that this small deficit in power at 
k ^ 0.2 h Mpc~^ is a numerical effect, most likely the 
transients in low-resolution simulations. 

How do the predicted non-linear power spectra depend 
on the amplitude of matter fluctuations? As the non- 
linear contributions to the power spectrum are given by 
the linear spectrum squared, a non-linear to linear ratio 
grows in proportion to cr| . In Fig2|we show how the non- 
linear contributions increase as one increases as from 0.7 
to 1.0. This figure may be useful when one compares our 
results with the previous work that uses different values 
of erg- 

5. DISCUSSION AND CONCLUSIONS 

The next-to-leading order correction to the matter 
power spectrum calculated analytically from 3rd-order 

APPENDIX 
CONVERGENCE TEST 

To test convergence of the power spectra derived from simulations and determine the valid range in wavenumber 
from each simulation box, we have run A^-body simulations with four different box sizes, Lbox = 512, 256, 128, and 
64 h~-^ Mpc, with the same number of particles, A^ — 256^. The initial redshifts are Zinitiai = 27, 34, 42, and 50 for 
512, 256, 128 , and 64 h^^ Mpc simulations, respectively. 

Figure lAll shows that simulations with a larger box size lack power on larger scales due to the lack of resolution, 
as expected, while they have better statistics on large scales than those with a smaller box size. This figure helps 
us to determine the valid range in wavenumber from each simulation box. We find that one can use 512, 256, 128, 



PT provides an almost exact description of the matter 
power spectrum in real space in the weakly non-linear 
regime, where A^(A:) < 0.4 (Fig.|2Jl. The most important 
implications of our results for the planned high-z galaxy 
surveys are that we can use PT to calculate (a) non- 
linearity in the baryonic acoustic oscillations (Fig. OJ, 
which should reduce systematics in constraining dark 
energy properties, and (b) the matter power spectrum 
up to much higher k than that was accessible before, 
which should vastly increase our ability to measure 
the shape of the primo rdial power spectru m as well as 
the mass of neutrinos l|Takada et al. 20(M) . Of course, 
these surveys measure the galaxy power spectrum in 
redshift space; thus, the future work should include 
PT calculations of non-linearity in (a) redshift space 
distortion (Scoccimarro 2004|). and (b) halo biasing 
(jFrv fc Gaztanaea 1993t iHeavens et al. 19 981. as weU 
as an extensive comparison with numerical simula- 
tions. PT also allows one to calculate the higher-order 
statistics such as the bispectrum, which has been 
shown to be a powerful tool to check for systemat- 
ics in our understanding of non-linear galaxy bias 
l|Mata,rreseeta,]. 19971 iVerde et a,]. 199^. We should 
therefore "reload" cosmological perturbation theory and 
make a serious assessment of its validity in light of the 
planned high-z galaxy surveys constraining properties 
of dark energy, infiation, and neutrinos. 

We would like to thank D. Ryu for letting us use his 
TVD code, and K. Gebhardt, Y. Suto and M. Takada for 
comments. D.J. would like to thank K. Ahn for his help 
on the TVD code. E.K. acknowledges support from an 
Alfred P. Sloan Fellowship. The simulations were carried 
out at the Texas Advanced Computing Center (TACC). 



and 64 ft, ^ Mpc simulations to calculate reliable estimates of the power spectrum in fc < 0.24 h Mpc 
0.5 h Mpc~\ 0.5 < fc < 1.4 h Mpc"\ and 1.4 < fc < 5 /i Mpc"\ respectively 



0.24 < fc < 
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Fig. A1. — Convergence test. {Left) Fractional differences between the power spectra from Af-body simulations in L^ox = 512, 256, 
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